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ORIENTATION EFFECTS ON THE INNER REGION OF DUSTY TORUS OF ACTIVE GALACTIC NUCLEI 
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ABSTRACT 

A sublimation process governs the innermost region of the dusty torus of active galactic nuclei. 
However, the observed inner radius of the torus is systematically smaller than the expected radius 
by a factor of ~ 1/3. We show that the anisotropy of the emission from accretion disks resolves this 
conflict naturally and quantitatively. An accretion disk emits lesser radiation in the direction closer 
to its equatorial plane (i.e., to the torus). We find that the anisotropy makes the torus inner region 
closer to the central black hole and concave. Moreover, the innermost edge of the torus may connect 
with the outermost edge of the disk continuously. Considering the anisotropic emission of each clump 
in the torus, we calculate the near-infrared flux variation in response to a UV flash. For an observer 
at the polar angle 8 ^ s = 25°, the centroid of the time delay is found to be 37% of the delay expected 
in the case of isotropic illumination, which explains the observed systematic deviation. 
Subject headings: accretion, accretion disks — dust, extinction — galaxies: active — galaxies: struc- 
ture — infrared: galaxies 



1. INTRODUCTION 

Active galactic nuclei (AGNs) are powered by gas ac- 
cretion onto supermassive black holes (BHs) at the center 
of each galaxy. The central BH and the accretion disk 
are hidden from our line-of-sight, in some geometry, by 
an optically- and geometrically-thick dusty torus (An- 
tonucci & Miller 1985; Miller & Goodrich 1990). Since 
the torus potentially plays a role of a gas reservoir for the 
accretion disk, its nature has long been studied (Pier & 
Krolik 1992; Efstathiou & Rowan-Robinson 1995; Beck- 
ert & Duschl 2004; Schartmann et al. 2008). 

A large geometrical thickness of the torus is inferred 
from the detection of high optical polarization, the sharp- 
edged ionization cone, the near and mid-infrared (NIR 
& MIR) to UV luminosity ratio and the number ratio 
between type-1 and type-2 AGNs in the unified AGN 
scheme (Antonucci 1993; Wilson & Tsvetanov 1994; 
Lawrence 1991). To sustain the geometrical thickness 
of the torus, the vertical velocity dispersion must be as 
large as ~ 100km s _1 . However, it is impossible to ex- 
plain such a large velocity dispersion by thermal velocity, 
since dust grains cannot survive over a sublimation tem- 
perature T su b ~ 1500K (Barvainis 1987; Laor & Drainc 
1993). Therefore, Krolik & Begelman (1988) concluded 
that numerous dusty clumps, rather than a smooth mix- 
ture of gas and dust, with a temperature < 1500 K con- 
stitute the torus with a large clump-to-clump velocity 
dispersion (i.e., dusty clumpy torus). Various models for 
the clumpy torus have been investigated (Nenkova et al. 
2002, 2008; Wada & Norman 2002; Dullemond k van 
Bemmel 2005; Honig et al. 2006). 

The dusty and clumpy torus absorbs optical/UV ra- 
diation from the central accretion disk, and re-emits as 
IR radiation (Telesco et al. 1984; Radovich et al. 1999). 
The energy balance (of each clump) between heating by 
incident radiation flux and cooling via blackbody radi- 
ation indicates that clumps located closer to the cen- 
tral BH have higher temperature. Clumps at the inner- 
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most region of the torus have the highest temperature, 
with a temperature of T su t, at the irradiated surface, and 
emit NIR radiation as "3/xm bump" (Kobayashi et al. 
1993). Barvainis (1987) derived the innermost radius of 
the torus (dust sublimation radius, denoted as i? su b,o m 
this study): 
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where Ltjv & n d a are UV luminosity and the size of dust 
grains, respectively. If this estimation is correct, the NIR 
radiation flux varies with a time lag of some months be- 
hind the optical/UV flux variations. 

Photometric monitoring observations of type-1 AGNs 
revealed that NIR emission from AGNs indeed lags be- 
hind optical variation by an order of a month (Clavel et 
al. 1989; Glass 2004; Minezaki et al. 2004). Moreover, 
the luminosity dependency of the time lag also coincides 
with the theoretical prediction as cx L^'v (Suganuma et 
al. 2006; Gaskell, Klimek & Nazarova 2007 using the 
data presented by Glass 2004). These agreements be- 
tween the observational results and the theoretical con- 
siderations support the idea that the innermost radius of 
the torus is controlled by the dust sublimation caused by 
irradiation from the accretion disk. 

Despite the successful agreement of luminosity depen- 
dency of the NIR-to-optical time lag, however, a conflict 
appeared concerning the normalization of the time lag. 
Namely, the time lag measured and collected by Sug- 
anuma et al. (2006) is systematically smaller than the 
lag predicted from eq. [T]by a factor of ~ 1/3 (Kishimoto 
et al. 2007; Nenkova et al. 2008). 

In this Letter, we notice that the estimation of dust 
sublimation radius (eq. Q]) presumes isotropic emission 
from the accretion disk. Yet, emission from an optically- 
thick disk is, in principle, anisotropic, which is a fact 
missing in previous theoretical considerations. There is 
a systematic difference between the inclination angle at 
which we observe the disk in type-1 AGNs and the an- 
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gle at which an aligned torus observes the disk. We 
show that the effects of the anisotropic illumination flux 
from the accretion disk resolve the systematic discrep- 
ancy between the observations and the theory regarding 
the NIR-to-optical time lag. In the next section, the 
anisotropy of the disk emission is introduced. Then, the 
inner structure of the torus is examined within the frame- 
work of clumpy torus models in §3. In §4, we derive the 
transfer function for NIR emission of the torus in re- 
sponse to optical/UV flux variation in the disk. Finally, 
we make a summary and discussion of this study in §5. 

2. ANISOTROPIC EMISSION OF ACCRETION DISK 

Now, let us consider the emission from an optically- 
thick, plain slab (i.e., disk). Radiation flux (F) from a 
unit surface area of the disk toward a unit solid angle 
at the polar angle of 9 decreases with an increasing 9 as 
follows: 



F oc cos 9 (1 + 2cos( 



(2) 



Here, the first term represents the change in the pro- 
jected surface area, while the latter represents the limb 
darkening effect for plasma, whose opacity is dominated 
by electron scattering over absorption (Netzer 1987 and 
references therein). In other words, an accretion disk 
emits lesser radiation in the direction closer to its equa- 
torial plane (i.e., larger 9; Laor & Netzer 1989; Sun & 
Malkan 1989). Consequently, the assumption of isotropic 
emission from accretion disks (e.g., eq. [T]) obviously over- 
estimates the radiation flux toward the torus. Therefore, 
it overestimates the inner radius of the torus. 

We note that this effect works even if the disk is in- 
finitesimally thin. The fact that the disk has a nonzero 
thickness brings about another anisotropy of illumina- 
tion flux, such that the torus is not illuminated below 
the disk height (near the equatorial plane) at 9 larger 
than a critical angle # m ax- Throughout this Letter, we 
adopt a thin disk (with an aspect ratio of ~ 0.01), like 
the standard accretion disk model (Shakura & Sunyaev 
1973). The effects of varying disk thickness as a function 
of accretion rates (Abramowicz et al. 1988; Fukue 2000; 
Kawaguchi 2003) will be investigated in a forthcoming 
paper. 

The effect of anisotropic emission has already been dis- 
cussed and referred to as the orientation effect in the con- 
text of emission lines, such as the Baldwin effect (Netzer 
1985; Francis 1993; Bottorff et al. 1997). However, con- 
sequences of this effect upon the torus have not been ex- 
amined so far. Hereafter, we will explore how this effect 
influences the torus structure. 

3. INNER STRUCTURE OF DUSTY TORUS 

We determine the inner edge of the torus so that the 
temperature of a clump (at the irradiated surface) equals 
to the sublimation temperature at the edge. As men- 
tioned above, radiation flux from the accretion disk varies 
with the polar angle 9. Thus, the sublimation radius of 
the torus is also a function of 9, which is indicated by 
Rsub(9) in this study (Figure HJ. Namely, R su b(9) is the 
distance between the torus edge and the central BH for 
various 9. To avoid confusion, we express the sublimation 
radius estimated under the isotropic emission assumption 
(i.e., eq. [1]), in contrast, as i? su b,o- 
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Fig. 1. — Schematic view of the inner structure of the torus. On 
the left side of the rotational axis, the definitions of £> m i n (the torus 
thickness) and 6> max (the disk thickness) are indicated. 
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Fig. 2. — Calculated geometry of the innermost region of the 
dusty torus, where r and z are the cylindrical coordinates. The 
central BH and the accretion disk in the 2 = plane are located 
at the origin (not drawn) and are observed with the polar angle 
Sobs = 25° . We assume alignment between the disk and the torus. 
In the case of an isotropic illumination, the torus edge has been 
supposed to stand at r = i? su b o ■ Thick solid line indicates the edge 
of the torus, with the opening angle of the torus (0 m itx) assumed 
to be 45°. On the right-hand side of the line, clumps contain dust. 
Dashed line represents R BUO (0) (eq. [3|. The torus inner edge is 
located closer to the central BH than in the previous estimations 
(eq. [TJ, and concave/hollow. The innermost region of the torus 
may connect with the outermost radius of the disk. 

Adopting the anisotropic illumination given in eq. [2j 
we obtain 
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Here, # bs is the polar angle toward the observer seen 
from the central accretion disk, and o bs = 25° is as- 
sumed throughout this study. Although various grain 
sizes result in the sublimation process occurring over a 
transition zone rather than a single distance (Nenkova et 
al. 2008), we employ a sharp boundary for simplicity. 
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Figure [5] shows the calculated structure of the inner- 
most region of the torus. The central BH and the accre- 
tion disk on the z = plane are located at the origin of 
the coordinate axes (not drawn). If the torus is indeed 
a reservoir of gas for the disk, angular momenta of the 
inf ailing gas will align these axes. Thus, we assume align- 
ment between the disk and the torus. In the case of an 
isotropic emission from the disk, the torus edge has been 
supposed to stand at r = i? su b,o- Thick solid line indi- 
cates the edge of the torus, with the opening angle of the 
torus (# m in) assumed to be 45°, and the maximum 9 of 
the torus (9 max ) set to 89° (i.e., thin disk approximation 
mentioned in the previous section). On the right-hand 
side of the line, the temperatures of clumps are below 
the dust sublimation temperature. Dashed line repre- 
sents i? S ub {0) (eq. |3]) and is drawn as a guide for various 
$min- It turns out that (i) the torus inner edge is located 
closer to the central BH than suggested by previous es- 
timations (eq. [1]) and that (ii) the structure of the edge 
is concave /hollow. These are the results of weaker illu- 
mination flux toward larger 9. Namely, dust can survive 
closer to the central BH at larger 9. 

Moreover, (hi) i? su b(#) decreases down to 0.1xi? su b,o at 
9 = 88.5°. This radius coincides with the outermost ra- 
dius of AGN accretion disks, which is determined by the 
onset of radial self-gravity of the disk (see Kawaguchi, 
Pierens & Hure 2004). Little is known about the re- 
gion outside the outermost radius of the disk (e.g., Collin 
2001), such as what happens in the 1-dex gap between 
the disk outermost radius and the torus innermost ra- 
dius (i?sub,o)- Our result indicates that there is no gap 
between the torus and the disk (Emmering, Blandford & 
Shlosman 1992; Elitzur & Shlosman 2006). 

4. TRANSFER FUNCTION 

Although interferometric NIR observations have been 
made for nearby Seyfert galaxies, only the visibility data, 
rather than images, are achieved because of the poor u- 
v coverage (e.g., Kishimoto et al. 2009). Even with the 
next generation telescope such as Thirty Meter Telescope 
having 0.01" spatial resolution at NIR, the innermost 
radius of the torus discussed above cannot be spatially 
resolved. Thus, observations of time variability will con- 
tinue to be powerful tools to explore the innermost struc- 
ture of the torus in the coming decade. In this section, 
we calculate the time variation of NIR emission in re- 
sponse to a (5-function like variation of the irradiation 
optical/UV flux (the transfer function). 

Transfer functions *&(t) for various geometries of the 
re-emitting region have long been studied, mainly in the 
context of the broad emission line region of AGNs (e.g., 
Netzer 1990). For instance, a ring produces a double- 
horned ^(t), while a thin (or thick) shell gives a rect- 
angle (or trapezoid). Time variation of the reprocessed 
radiation (NIR in this study) is a convolution of the il- 
lumination flux variation with ^(t). The measured time 
lag corresponds to the centroid of ^(t). Below, we calcu- 
late ^{t) and its centroid in order to compare them with 
the observational results. 

It turned out that it took ~ 1 year for the inner re- 
gion of the dusty torus to adjust to the varying illumi- 
nation flux (Koshida et al. 2009; Pott et al. 2010). We 
hereafter regard the inner structure of the torus as time- 
independent in the timescale of NIR-to-optical time lag 



(~months). 

4.1. Calculation Method 

To calculate ^(t) for the clumpy torus, we discuss the 
following three items; (a) the optical path, (b) NIR emis- 
sivity of the torus inner region and (c) anisotropic emis- 
sion of each clump. In this Letter, we ignore the re- 
sponse from the torus edge at 9 > -|, assuming torus self- 
occultation (absorption of NIR emission from a clump by 
other clumps in the line of sight). Although we do not 
go further, this assumption can be tested in principle via 
the profile of broad emission lines and its time variation 
(e.g., Peterson 2001). 

Firstly, (a) the optical path difference is written as 

-Rsub(#) [1 - {cos 9 ohs cos 9 + sin 9 ohs sin 9 cos <j)}] , (4) 

where 4> is the azimuthal angle and defined so that <f> — 
for the observer (Figure [T]). The concave shape of the 
inner region of the torus reduces the optical path differ- 
ence. Clumps at a distance slightly larger than i? su b(0) 
with a slightly lower temperature than the sublimation 
temperature will somehow emit NIR radiation. More- 
over, clumps at a distance slightly smaller than i? su b(#), 
in the middle of the dust destruction process but with 
some of their dust grains still surviving, will also emit 
NIR radiation. These effects will smear out the resultant 
NIR response ^f(t). Since these effects are unlikely to 
change the centroid of ^f(t) drastically, we here ignore 
them and consider only optical paths that hit the inner 
edge of the torus. 

Next, (b) we discuss the emissivity of the torus in- 
ner region as a function of 9. Since i? su b(#) varies 
with 9, the distance dependency of the NIR emissiv- 
ity must be determined. Following the calculations of 
three-dimensional radiative transfer (e.g., Honig et al. 
2006; Schartmann et al. 2008), we assume that (b-1) the 
clump size increases and (b-2) the clump number den- 
sity decreases when the clump-to-BH distance increases, 
as follows. 

Closer to the BH, the tidal force is stronger, thus only 
smaller clumps are likely to survive there. Therefore, 
we assume that the clump size is in proportion to the 
distance between the clump and the central BH. Then, 
the apparent angular size (observed from the BH) of the 
clumps at the torus inner edge, which are emitting NIR 
radiation, is independent of 9. (Farther clumps receiving 
direct emission from the disk as well are not of interested 
in this study, because they emit MIR rather than NIR.) 

Little is known about the number density of clumps. 
We assume that the number of clumps at the torus edge 
per unit solid angle (g?£!) is independent of 9. Namely, 
the same number of NIR-emitting clumps is contained 
per d£l for various 9. To summarize the two factors, the 
emissivity of NIR flux per dQ is assumed to be propor- 
tional to i? su b(60 2 . 

Lastly, (c) the anisotropy of NIR emission from each 
clump is considered. To calculate NIR flux, we must care 
about how extent the illuminated surface of a clump is 
seen by the observer (Nenkova et al. 2002). Honig et al. 
(2006) mentioned that this effect resembles the phases 
of the moon. Let us suppose that an observer looks at a 
clump with an angle £, where £ = means a face-on view 
of the illuminated surface (the angle denoted as <j> in the 
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Fig. 3. — Resultant transfer function for the torus we consider 
(solid line; NIR response to a 5-function like variation of illumi- 
nation flux). Dashed line is the transfer function for a ring with 
radius i? su b,o an d $obs = 25°. In the case of the anisotropic disk 
emission, the torus responds approximately three times faster than 
the ring does (the centroid of the transfer function is 0.37ii su t, q/c). 



paper by Honig et al.). In the case of the moon phase, 
the anisotropy coefficient (the fraction of the illuminated 
surface seen by the observer) is (1 + cos£)/2. However, 
incident photons can somehow penetrate gaseous clumps, 
heating up the limb of each clump. Then, the anisotropy 
is weakened compared with the above estimation for the 
solid (moon). We adopt the following anisotropic coeffi- 
cient for the waning effect, 



mm 
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1 + cos £ 



0.1 



(5) 



This coefficient is chosen so as to reproduce the Monte 
Carlo calculations by Honig et al. (2006) for a single 
clump observed from three different £. We also try other 
descriptions for the anisotropy under a condition that 
each description reproduces the Monte Carlo results, ob- 
taining little difference in the resultant $f(t). 

4.2. Result 

Figure [3] presents the resultant transfer function of the 
torus for 9 b s = 25° (solid line) , calculated by integrating 
9 (from # m i n to ma ,x) and <f> (0 to it). The portion at 
cj) ~ produces the left horn, which is a faster (due to a 
shorter optical path difference, eq. U) and weaker (due 
to the waning effect, §4.1) response than the opposite 
area at <fr ~ 7r (right horn). Dashed line is the transfer 
function for a ring with radius i? S ub,o observed from 9 Q b s 
= 25°and is drawn for an example of transfer functions in 
the case of isotropic illumination flux. The ring responds 
to the illumination flux at around R S ub,o/c (where c is 
the speed of the light) after the (5-function like variation 
in the illumination flux at zero time lag. (It is just a 
coincidence that the ring starts to respond exactly when 
the torus ends with the current parameter sets.) 

On the other hand, the torus we consider shows a faster 
response and narrower width of the transfer function. 
These results are due to the facts that the torus inner 
region we consider is (i) closer to the central BH than 
-Rsub,o and (ii) concave. The centroid of the time delay 
of the transfer function is found to be 0.37i? sub ,o/c, which 



solves the puzzle of the systematic deviation of a factor 
of ~ 1/3 between the observational results and the theo- 
retical considerations (under the assumption of isotropic 
emission) . 

5. SUMMARY AND DISCUSSION 

The dusty clumpy torus surrounds the accretion disk 
and BH in AGNs. Various observations have revealed 
that the dust sublimation process governs the innermost 
region of the torus. However, there was a systematic de- 
viation between the observational results and the theory 
regarding the inner radius of the torus. 

In this study, we have shown that the anisotropy of the 
emission from accretion disks resolves this conflict natu- 
rally and quantitatively. Namely, the angle at which we 
observe the disk in type-1 AGNs is systematically closer 
to a pole-on view than the angle at which an aligned 
torus observes the disk. An accretion disk emits lesser 
radiation in the direction closer to its equatorial plane 
(larger 9). We have found that the anisotropy makes 
the torus inner region closer to the central BH and con- 
cave. Furthermore, the innermost edge of the torus may 
connect with the outermost edge of the accretion disk 
continuously. 

Considering the anisotropic emission of each clump, 
we have calculated the NIR flux variation of the torus in 
response to a UV flash. A rapid response to the illumi- 
nation flux is realized, due to the small inner radius and 
concave shape of the torus. For an observer at o bs = 25°, 
the centroid of the time delay is found to be 37% of the 
delay expected in the case of isotropic illumination. 

Other than the orientation effect examined in this 
study, large grain size (see eq. [1} and/or extinction be- 
tween the torus and the accretion disk (i.e., in the broad 
emission line region) are possible concepts to reduce the 
inner radius of the torus (Laor & Draine 1993; Maiolino 
et al. 2001; Gaskell et al. 2007). 

Anisotropy is an inevitable feature of the emission from 
optically-thick disks. Even if a BH object (seen in nearly 
face-on geometry) looks very luminous around the Ed- 
dington luminosity, the flux in the directions close to the 
equatorial plane is not necessarily huge to prevent gas 
infall. This effect (along with the disk self-occultation) 
enables gas supply to luminous BH objects. Actually, the 
observed data do not support the concept of Eddington- 
limited accretion (Collin & Kawaguchi 2004). 

Future studies will examine misalignments between the 
torus and the disk, and different o bs (for type 1.5, 1.8 
etc objects), 9 mm (for objects with a thick/thin torus) 
and 9 max (for thick disks with super-Eddington accretion 
rates). For instance, some Narrow-Line Seyfert 1 galax- 
ies and Narrow-Line QSOs prohibit weak NIR emission 
(Rodrfguez-Ardila & Mazzalay 2006; Kawaguchi et al. 
2004; Jiang et al. 2010; however, see Hao et al. 2010). 
Small # m ax due to the disk self-occultation with super- 
Eddington accretion rates (Fukue 2000) can be a reason 
for the weakness. 
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